Dynamical mean-field theory of electron-phonon interaction in correlated electron 
materials: general results and application to doped Mott insulators 
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The dynamical mean-field method is used to formulate a computationally tractable theory of 
electron-phonon interactions in systems with arbitrary local electron-electron interactions in the 
physically relevant adiabatic limit of phonon frequency small compared to electron bandwidth or 
interaction scale. As applications, the phonon contribution to the effective mass of a carrier in a 
lightly doped Mott insulator is determined and the phase separation boundary is discussed. 
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'Correlated electron' materials such as doped Mott in- 
sulators ffl, high-temperature superconductors || , heavy 
fermions |J , and 'half- metallic' oxides Q are central to 
present-day condensed matter physics because they ex- 
hibit many properties that seem incompatible with pre- 
dictions of the standard paradigm [local density approx- 
imation (LDA) band theory plus Boltzmann transport] 
of metal physics. A crucial role is played by electron- 
electron interaction effects beyond the scope of the LDA, 
and theoretical attention has focused on these effects ne- 
glecting other physics. However, the electron-phonon in- 
teraction is present in all real materials, and it is therefore 
important to understand its effects in correlated systems. 
The relevance of electron-lattice effects in the CMR man- 
ganites is by now well established (|J|. Lanzara and 
co-workers presented photoemission data which, they ar- 
gue, indicate that electron-phonon effects play an impor- 
tant role in high- temperature superconductivity Q. On 
the other hand, the apparent absence of electron-phonon 
effects in the resistivity of high-T c materials is a long- 
standing mystery 0. 

There is to date no systematic theoretical extension to 
correlated materials of the successful Migdal-Eliashberg 
(ME) theory S, which describes electron-phonon inter- 
actions in weakly correlated materials. Electron-phonon 
interactions in one-dimensional systems have been stud- 
ied by renormalization group methods |}| , and recent im- 
provements in numerical techniques have for example al- 
lowed the spin-phonon dynamics of insulating quasi-one- 
dimensional Peierls systems to be determined in consid- 
erable detail pX([ . Concerning higher-dimensional, metal- 
lic systems, Kim and co-workers |ll[] argued that near a 
Mott transition a decrease in electron-phonon coupling 
was compensated by an increase in carrier mass, leading 
to electron-phonon effects unrenormalized by proximity 
to the Mott transition. Our results, to be presented be- 
low, disagree with this conclusion. Several authors have 
used direct numerical simulation [e.g., quantum Monte 
Carlo (QMC) techniques] of models of electrons inter- 
acting both with each other and with lattice vibrations 



]l2| [l4|. However, the present limitations of memory size 
and algorithms have restricted these works mainly to the 
calculation of static properties, especially phase bound- 
aries and transition temperatures and to the 'antiadia- 
batic' limit of phonon frequency comparable to electronic 
energy scale. 

The development of the dynamical mean-field (DMF) 
method Jl5[ has opened an important avenue for 
progress, by showing that if (as occurs for the electron- 
phonon interaction) the momentum dependence of the 
electron self-energy is negligible then a good approxi- 
mation to the correlation physics can be obtained from 
the solution of a numerically tractable quantum impurity 
problem plus a self-consistency condition. Unfortunately 
the straightforward inclusion of the electron-phonon cou- 
pling in the DMF formalism is difficult because the 
mismatch between the typical phonon frequency scale 
uiq < 0.1 eV and electron energy scale t > 1 cV ren- 
ders conventional numerical approaches to the impurity 
problem prohibitively expensive, except in the 'antiadia- 
batic' limit Jl4| of relevance to rather few materials. In 
this paper we present a practical implementation of an 
adiabatic expansion of the DMF formalism and show how 
it may be used to determine the phonon contribution to 
electronic properties of correlated systems. The work re- 
ported here builds on previous papers, which introduced 
the adiabatic expansion JLsf and applied it to models in- 
volving only electron-lattice interactions p7| , |lq] . As an 
application we answer the long-standing question of the 
electron-phonon contribution to the carrier self-energy in 
a lightly doped Mott insulator. Our methods may easily 
be extended, e.g., to heavy fermion materials and may 
be combined with recent extensions of the DMF method 

We consider the single-site DMF approximation to a 
general Hamiltonian of the form H = -ffband + ^ci-ci + 
H p h + -ffci-ph, where -ffband an d -Hei-ei describe electrons 
moving in some band structure and interacting via some 
local interaction. We model the phonons as dispersion- 
less quantum oscillators with instantaneous displacement 
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q a , mass M a , and spring constant K a : 



(1) 



Dispersion may be approximately included by treating 
M, K as Brillouin-zone averages of phonon dispersions 
or via the more sophisticated techniques of Refs. [|TJ and 
\22\. Anharmonic terms such as those considered in Ref. 
[ 23 1 can be easily added and will be seen to be generated. 

The electron-phonon coupling is again taken to be lo- 
cal. It is cumbersome to write in general. We illus- 
trate the issues via the physically relevant example of 
the pseudocubic manganese perovskites. Here the rele- 
vant electrons are e 9 -symmetry Mn-0 hybrid states. The 
generalized phonon field q a encompasses (i) a 'breathing' 
mode (symmetric distortion of the Mn-Oe octahedron) 
coupling via a constant gB to the total on-site charge 
density and (ii) a one-parameter family of Jahn- Teller 
modes (even-parity, volume-preserving distortions of the 
Mn-Og octahedron) coupling via a matrix element gjT 
to appropriate differences of occupancy between different 
orbitals (for a more detailed discussion of the couplings 
and notation see, e.g, Ref. j|]). These two distortions will 
be labeled by a scalar coordinate x and a two-component 
vector Q = (Q x ,0,Q z ), respectively. The coupling term 
reads (repeated indices are summed) 
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/}gBXj(nj -n)+ gjrQi 



-iafi ' 
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(2) 



where r = (t x , t v ,t z ) are Pauli matrices acting in orbital 



n.a C iaa C iaa 



is the local electron density, 



space, rii 

and (xi) = is defined to be the equilibrium phonon 
state for a uniform distribution of electrons. 

Within the single-site DMF approximation |D| the 
properties of H may be obtained from the solu- 
tion of an impurity model specified by the action 
S[c,c,q,a] = S [q] + Si[c,c,q,a], with S [q] = 
^I^T)Y Jk dk{K a +M a ul)q a _ k and Si fee,?, a] = 

S e l_el[c, c] + S c i_ p h[c, C, q] - J2 a „ 

Canaan a n . Here the 



terms S c i-ci and S i_ p h are obtained in the usual way 
from the interaction terms listed above, and S is a func- 
tional of a mean-field function a, which may depend on 
spin and orbital indices and expresses the effect of the rest 
of the lattice upon the single site. It is fixed by equating 
the impurity Green function Gi mp [a]„ = 5hxZ[a]/5a n to 
the local Green function Gi oc [a]n = / d 3 p/(27r) 3 (iu) n + 
fi-H[a] n - -Hband) -1 with £[o]„ = c - G imp [a]" 1 . The 
fields c, c and q represent local electronic and phonon 
degress of freedom and have been Fourier transformed 
according to q(r) — J2k ex P(~ * w fe r )% e ^ c - We will index 
bosonic Matsubara frequencies ujk — 2knT by integers k 
and fermionic Matsubara frequencies lu h = (2n + 1)ttT 
by integers n, i.e., q k = q{iu k ) and a"* = o" 1 ^^). 

To analyze the effect of phonons on electronic physics 
p6| , p7[ we integrate out the electron fields and work with 



an effective action S[q, a] — So[q] + Si [q, a], which we ex- 
pand about the values q that extremize S. We introduce 
an electronic bandwidth or interaction scale t and define 
the parameters 



7 



{K/Mf/ 2 



t 



t ' 



Crucial objects in the expansion are the vertices 

1 5Si[q,a] 



rjv[o]/s 1 ,... ) fej, 



(N-l)\ 5q a k \---8qZ 



(3) 



(4) 



which are connected correlation functions of the 
electrons-only theory. The adiabatic parameter 7 con- 
trols the expansion because each phonon loop involves a 
sum over frequencies of order loq ~ jt while the vertices 
r vary with frequency on the scale t. If terms of order 
7 3 / 2 and higher are neglected, the IV^ may be dropped 
and 1^3,4 may be approximated by their static value com- 
puted using a at 7 = 0. With suitably rescaled fields and 
frequencies (for details see Ref. [H) we may write 



S[o, a] = )- ^ 1k D k V 



A a /a i/ayi/a 

3 r 3H 2^ ^q^q-^-k 



A 2 7 T 



r 4 [a] ^2 Qkiqk 2 <ik 3 q-k 1 -k2-k 3 , (5) 



ki,k 2 ,k 3 



where D k is the phonon propagator, equal for scalar 
phonons to 1/(1 + uj^/uJq — \T2[a]k.-k), renormalized 
by electron-electron interaction effects contained in IV 
Note that 0(7) corrections to a, as well as the leading 
frequency dependence, must be included in IV 

We used the Hirsch-Fye |24j algorithms distributed 
with Ref. [lj| and routines to compute the connected 
correlation functions for a single-band Hubbard model. 
Computations were performed on a Sun workstation; the 
longest computation T2 at larger U took ~ 40 minutes. 
Errors (estimated by comparison to exact solution and 
scaling with system sizes) were less than 2% for param- 
eters shown. Fig. |l| shows that near half filling and 
for scalar phonons the Hubbard-JJ strongly decreases T2, 
thereby suppressing the polaronic instability occurring 
at = 1. Although the polaronic instability is sup- 
pressed a phase separation instability is found for A > 1 
and n near 1; details will be presented elsewhere. For 
comparison and to demonstrate the power of the method 
we also present in the right-hand panel results for a model 
of 'colossal magnetoresistance' manganites in which or- 
bitally degenerate electrons are coupled to Jahn- Teller 
(JT) phonons via the second term in Eq. (|^) and also 
feel both a Hubbard-C/ and a 'double-exchange' interac- 
tion arising from a strong coupling to core spins (i.e., the 
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model defined in Ref. H but with a Hubbard-JJ added). 
The U is seen to enhance the effect of JT phonons. Fig. 
H shows that for the Hubbard model and T4 are even 
more rapidly suppressed by U so repulsive interactions 
effectively decouple electrons from scalar phonons. 
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FIG. 1. Left panel: quadratic vertex F2 for scalar phonons 
as a function of electron density n, for various values of elec- 
tron-electron interaction u := U/(2t) and /3 — 10. Right 
panel: T2 for Jahn- Teller phonons in the paramagnetic 
(i.e., spin-disordered) phase of the Hubbard-double-exchange 
model -ff c i_ c i = Un(n — 1) + JS ■ a with J — > oo. Numerical 
method: QMC with L — 32 time slices. 




FIG. 2. Left panel: cubic vertex T3 for single-band Hub- 
bard model at f3 = 10, with u and electron density n as shown. 
Right panel: quartic vertex T4 for same parameters. Numer- 
ical method: exact enumeration with L = 16. Note: f3 = 10 
is higher than the critical endpoint of the Mott transition 
occurring at u ~ 1.5. 

We now turn to the electron self-energy E ph due to 
interactions with phonons, assuming that the effective 
vertices r^v have been determined as discussed above. 
The calculation is simplified and made more physically 



transparent if the action is viewed as a functional of 
the full electron propagator G rather than as a func- 
tional of a. The key object is then the 'local Luttinger- 
Ward functional' 4>[G] defined 25 as the sum of all 
vacuum-to-vacuum skeleton diagrams and related to the 
self-energy via S = 5(f>[G]/6G. Within DMF theory, 
<j> may be found from the local thermodynamic poten- 



tial fijmp = — TlnZ (see Refs. |15 and [E6|) so may be 



computed by an adiabatic expansion. We find <p[G] — 
</> cl - cl [G]+0 ph [G]+O(7 3 / 2 ) with (f>P h [G} = lEfeln-Dfe 1 . 
Therefore 



(6) 



In general T 2 and thus Sr 2 /5G vary on the scale 
of t so S ph ~ 7, an unimportant correction to 
the bare frequency or electron-electron contribution. 
However, if 5T/SG is singular at low frequencies, 
as in Fermi liquids, then dH/duj may be of order 
unity. To investigate this we note that r2[G]fc,-fc = 
-2tTj2 n G n G n+k (1 + 2tT£„, k* n , k G nl G n , +k ), where 
the particle-hole reducible vertex A R is given in terms 
of the particle-hole irreducible vertex A 1 via the Dyson 
equation shown in Fig. |§|. 
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3. Diagrammatic Dyson equation connecting the re- 
ducible and irreducible particle- hole vertices A R and A , re- 
spectively. 



The standard analysis implies that A 1 is a smooth 
function of its arguments, so the required singular be- 
havior can only occur if we differentiate on one of the 
explicit G factors. Further, to 0(7) we may neglect the 
frequency dependence of the vertices leading to 



E ph = XtTj2DkG n+k A 2 + S ph 



rcg 



(7) 



with S ph rcg varying with ui n on the scale of t and 

A = l + 2fT^A« o G£ + 0(7). (8) 

n 

If the ground state is a Fermi liquid then G{iv) = 
— iir sign(i') p(fi) + G mc , where G lnc is nonsingular at 
small frequencies [E7|. The leading contribution in 7 to 
m*/m is then m*/m\ph = 1 + Xtp(fi)A 2 + 0(7), where 
A = A/(1-Ar 2 ). 

A is a vertex function of the underlying electron- 
electron theory, determined in general by solving the ver- 
tex equation shown in Fig. ^. However, to the order 
to which we work it is simpler to note that in a Fermi 
liquid the leading low-frequency behavior of the density- 
density correlation function is x(w) = x(0) + iA\u>\ with 
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A = A 2 d X °/duj and xl = ~ 2tT T, n G n G n+k is the func- 
tion obtained by convolving two exact Green functions 
with no vertex corrections. Thus 
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,2 _ dxk/d(iuj k )\k=o 



A' = 



dx k /d(iuj k )\ k=0 



(9) 



The right hand side of Eq. @ can be calcu- 
lated numerically by QMC in the time domain be- 
cause the leading long-time behavior comes from the 
first frequency derivative and is x(0 <C r <C /3) = 
(Tr/P)dx k /d(iu k )\ k =o/s™ 2 (irT/{3). For /3 > 4 we find 
x(t)/x°(t) is very flat in the center region of the inter- 
val (0, (3) (see inset to Fig. ||). The rapid decrease of A 
again shows that the electron-phonon mass enhancement 
and scattering rate are 'turned off' near half filling and 
for strong correlations, in disagreement with the authors 
of Ref. (n[] who found a weaker decrease in the matrix 
element, which was moreover cancelled by an increase in 
effective mass. 
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FIG. 4. Static particle-hole scattering renormalization fac- 
tor as function of u := U/(2t). Numerical results for (3 = 16 
and n — 0. Circles: data as described in text. Dashed line: 
analytic weak-coupling approximation. Inset: imaginary time 
dependence of \ showing accuracy of approximation. 

We summarize the main results of this paper. We have 
studied a general electron-phonon model in a new ap- 
proach that combines the adiabatic expansion in 7 < 1 
of the conventional Migdal-Eliashberg theory with a dy- 
namical mean- field treatment of electron correlations. 
Our main finding is an effective phonon action correct to 
C(7 3 / 2 ) with static coefficients T N (N — 2, 3, 4) and ver- 
tex A, which are easy to compute numerically for various 
models. For the Holstein-Hubbard model we have shown 
that the electron-phonon coupling is strongly suppressed 
by Hubbard-t/ effects in systems near Mott transitions. 
Extension of our results to other correlated systems such 
as heavy fermions would be worthwhile. 

We thank S. Blawid for useful discussions, H. Monien 
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